Elastic and inelastic collisions of 6 Li atoms in magnetic and optical traps 
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We use a full coupled channels method to calculate collisional properties of magnetically or 
optically trapped ultracold 6 Li. The magnetic held dependence of the s-wave scattering lengths of 
several mixtures of hyperfine states are determined, as are the decay rates due to exchange collisions. 
In one case, we find Feshbach resonances at B = 0.08 T and B = 1.98 T. We show that the exact 
coupled channels calculation is well approximated over the entire range of magnetic fields by a simple 
analytical calculation. 
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The observation of Bose-Einstein condensation in 
atomic alkali gases has triggered an enormous in- 

terest in degenerate atomic gases. At present, one of the 
most important goals is to achieve quantum degeneracy 
in a fermionic gas. In the case of fermionic 6 Li, it has 
been shown theoretically that a BCS transition to a su- 
pcrfhnd state could be realized at a critical temperature 
on the order of temperatures obtained in the BEC ex- 
periments [il. This relatively high critical temperature 
is due to the fact that 6 Li has a very large and nega- 
tive triplet s-wave scattering length ar — — 2160ao ||, 
and that at sufficiently large magnetic fields a mixture of 
the upper two hyperfine states |6) and |5) is essentially 
electron-spin polarized 

The disadvantage of such a large triplet s-wave scat- 
tering length is that the exchange and dipolar relaxation 
rates for the gas are also anomalously large. Neverthe- 
less, to suppress this decay, one can apply a magnetic 
bias field. In Ref. Q], we used the distorted wave Born 
approximation (DWBA) to calculate the corresponding- 
decay rate constants for these decay processes and found 
that at large magnetic fields B > 10 T, the dipolar 
rates are dominant, but at smaller magnetic fields, the 
exchange rates greatly exceed those due to the dipolar 
interaction. However, the DWBA is only valid at mag- 
netic fields B > 0.1 T, and we were at that time unable 
to make predictions at lower, experimentally more con- 
vienent fields. 

The aim of the present paper is to provide useful in- 
formation on the s-wave scattering length and exchange 
decay rate constants at lower magnetic fields. In view 
of the ongoing experiments, we will concentrate on colli- 
sions involving the following antisymmetrized hyperfine 
states: |{65}>, |{64}), |{54}>, and |{21}>. The first three 
mixtures contain states that are low-field seeking at suf- 
ficiently high field, and therefore, can be confined in a 
magnetic trap. In contrast, the combination |{21}) can- 
not be magnetically trapped, but can be confined in a far 
off- resonance optical trap. We do not consider any other 
combination of high-field seeking states, since the |{21}) 
mixture cannot decay through collisions and is therefore 
most favorable experimentally. In addition, Van Abeelen 
et al. have already considered the |{62}) combination, 



which is low-field seeking at very weak magnetic fields 
B < 26 x 10~ 4 T, but does not have as large an s-wave 
scattering length as some other combinations ||. 

To obtain the magnetic field dependence of the s- 
wave scattering length and exchange decay rate constants 
in the various cases, we perform a full coupled chan- 
nels (CC) calculation M, using the most up-to-date sin- 
glet and triplet potentials Vo(r) and Vi(r) ||. Further- 
more, we show that the rate constants and the scattering 
lengths found using the CC calculation can be obtained 
analytically using a simple approximation that we call 
the asymptotic boundary condition (ABC) approxima- 
tion Jnj|. 

If the thermal energy is much smaller than the hyper- 
fine plus Zeeman energy, A Q/ g^ a /^/, gained in the tran- 
sition from the incoming state |{a/3}) to an outgoing 
state \{a'(3'}), the corresponding exchange rate constant 
is given by the zero temperature expression M 
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In this expression, fi is the reduced mass of the two 6 Li 
atoms, hk a /3 is the relative momentum of the incoming 
particles in state |{a/3}), and the matrix S is the scat- 
tering matrix of the multi-channel problem with angu- 
lar momentum quantum numbers I = and m = 0. 
The number of channels being coupled is determined by 
the fact that the central interaction V c {r) = Vo(r)P° + 
Vi(r)'P 1 cannot change the total nuclear plus electron- 
spin projection of the two-particle wave function along 
the magnetic field. So for example, the state |{65}) can 
only decay to |{61}). While the nondiagonal part of the 
unitary S-matrix gives the decay rates, the diagonal part 
of the S-matrix gives the s-wave scattering length via 

S{ a p},{ a p} — > e x P (-2jfc Q/3 a a/3 ) ~ 1 - 2ik aP a a p, , (2) 

when k a p — > 0. The goal of the CC calculation is to de- 
termine the S-matrix, but rather than giving a detailed 
description of the calculation, we refer to Ref. || and 
merely present the results here. 
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In Fig. |], the real part of the s-wave scattering length 
is plotted for the four cases of interest. 
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FIG. 1. Real part of zero-temperature s-wave scattering 
length for 1) 65, 2) 54, 3) 64, and 4) 21 collisions as a func- 
tion of the magnetic field. 

Note that, in the first three cases, the s-wave scat- 
tering lengths attain the large and negative value of 
a T = -2160a , only for B > 0.05 T. At lower mag- 
netic fields, the scattering lengths become 47ao at zero 
magnetic field for the first three cases, and zero in the 
latter case. This effect is explained in detail below. The 
scattering length 021 exhibits Feshbach resonances Jll|] at 
B ~ 0.08 T and B ~ 1.98 T. Physically, a Feshbach res- 
onance arises whenever the Zeeman energy of the incom- 
ing wave function |{21}), which is almost purely triplet, 
coincides with a bound-state energy of the singlet poten- 
tial Vo(r). Only the state |{21}) can exhibit Feshbach 
resonances at the magnetic fields of interest since it is 
the only combination that has a negative total energy 
(about —2/i e B). The precise location of these resonances 
is shifted slightly due to the hyperfine coupling, but their 
positions can be accurately computed since the binding- 
energies of the relevant singlet states, v — 38 and v = 37, 
are accurately calculated from the experimentally deter- 
mined singlet potential ||. 

The decay rate constants resulting from the CC cal- 
culation are, for all possible exchange processes, plotted 
in Fig. |^. The sharp minimum in the rate G^^n is a 
result of the fact that the matrix element ({41}|I /C |{54}) 
is zero if tan(9_ = l/y/2. However, the state |{54}) still 
remains coupled to |{41}) through the other elements in 
the (5 x 5) coupling matrix. As a result of this interfer- 
ence, the rate constant is never zero, and the location of 
the minimum is slightly shifted from the magnetic field 
at which tan = 1/ \[2. 

The second remark regarding Fig. |2| concerns the be- 
havior of the rate constants at large i?-fields where 
sin0± ~ 9± oc 1/B. The slope of each curve 
is determined by the magnetic field dependence of 



'B \ {{a' (3'}\V c \{a(3})\ . After some algebra, one can 
show that for curves 2, 4, 6 and 7, the relevant matrix 
element is equal to 6±(Vi — Vq)/2 = 9±V ex /2, which is a 
factor 1/2 smaller than in the cases 1 and 3. Therefore, 
the decay rates are a factor of 4 smaller than in the latter 
cases, but the slopes are identical. On the other hand, 
the matrix element ({21}|F C |{54}) = 39+6-V ex /2, has 
an additional factor of 1/B, and thus the slope of curve 
5 is a factor of 7/3 smaller than the others. 
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FIG. 2. Exchange decay rate constants for the processes 1) 
65 -» 61, 2) 54 -> 63, 3) 54 -> 52, 4) 54 -» 41, 5) 54 -> 21, 
6) 64 — > 62, 7) 64 — > 51 as a function of the magnetic field 
strength. 

A comparison of the rate constants shows that the 
I {64}) combination is slightly more favorable than the 
I {65}) combination. The total rate constant of the for- 
mer combination is only half the rate constant of the lat- 
ter at large magnetic fields. At B — 0.1 T, they amount 
to 1.49 x 10 -9 cm 3 /s and 2.74 x 10~ 9 cm 3 /s, respectively. 

We now show that the CC results can be accurately re- 
produced using a simple analytical calculation. We con- 
sider here only the 65 case, but calculation of the other 
combinations is analogous, and gives similar accuracy. In 
the ABC approximation, space is divided in two regions. 
In the interior region r < R, the singlet and triplet po- 
tentials dominate, while in the outer region r > R, the 
hyperfine energies prevail. The boundary R is chosen be- 
tween the points where the exchange potential V ex (r) is 
of the order of the hyperfine constant a^f, and the point 
where the triplet and singlet potentials themselves are 
of order a/j/. In our case this is between r — 28ao and 
r = 62ao, and below we will use as examples R = 40ao 
and R — 60ao. 

In the interior region the hyperfine splitting is ne- 
glected, so that the channel wave functions are just linear 
combinations of the singlet and triplet scattering wave 
functions, which we approximate by their asymptotic 
forms As smk a p(r — as) and At sm.k a p(r — ay), in the 
limit of small momentum k a p. The known values of the 
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scattering lengths are as = +45.5ao and ax — — 2160ao, 
respectively g. In the exterior region, we neglect the 
central interaction V c , so the channel wave function 
is a plane wave exp (—ik a pr) for the incoming channel 
\{a/3}), and #{ a '/3'},{ a( 3} exp {+ik a >p>r) for all outgoing 
channels |{a'/3 }). The wave numbers fc Q '/3' depend on 
the energy gained by the transition to the respective out- 
going channel, that is, k a > Si > = \Jk 2 aj3 + 2/iA ct/3 ^ Q / /3 //^ 2 . 

The yet unknown amplitudes As, At, and S \ a i piy ,{ a p} 
must be determined by imposing continuity and differen- 
tiability to the wave function in each channel at r = R. 
To do so, we need to find the precise linear combination 
of singlet and triplet wave function in the interior region. 
In terms of the basis IS" Ms; I Mi), where S = si + S2 
and I = ii + i2, are the total electron and nuclear spin 
of the two-atom system respectively, we find 

|{65}) = sin(9 + |0 0;2 2) +cos(9 + |l 1; 1 1), 
|{61}) =cos(9+|0 0;2 2) -sin0+|l 1; 1 1). 

Evidently, for the channel wave functions we have to take 
linear combinations with the same coefficients, so, for ex- 
ample the wave function of the |{65}) channel in the inte- 
rior region becomes ^65 i r ) — As sin 9+ sin k a p(r — as) + 
At cos #_|_ sin k a p(r — ax). It can be shown easily that 

are related with the S-matrix 



the coefficients S 
H according to 



{a'f3'},{af3} 
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For the case of the incoming |{65})-state, we are only 

interested in the constants >5{65},{65} ancl 5'{6i},{65}j an d 
therefore can eliminate the constants As and At from 
the 4 x 4 set of complex equations. The resulting 2x2 
set of complex linear equations for these amplitudes can 
be solved in orders of k§§, and we find to first order in 

&65 



s. 



{65},{65} 



= -l + 2ih 



65 
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ik^{R-a s ){R-a T ) 



with the CC curve, and the result of a degenerate inter- 
nal states (DIS) approximation ||. Fig. ||o shows the 
real part, while Fig. pp gives the imaginary part. 
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FIG. 3. The a) real and b) imaginary part of the s-wave 
scattering length a§§. The solid line is the CC value, which 
goes to the triplet scattering length for large magnetic fields. 
The dashed lines are 1) the DIS approximation, 2) the ABC 
approximation with R = 40ao, and 3) the ABC approxima- 
tion with R = 60ao. 



cos 2 9 + (R-a T ) + sin 2 6+ {R-a s ) 



where 



(4a) 



~ 2ik e5 cos 9+ sin 6 + (as - a T ) .. 
<!>{6l},{65} = ' ( 4b ) 



2/xA 65 _ (6 i/?i 2 , and D(k$) = 1 



ikgi [cos 2 0+(R — as) + sin 2 # + (i? — ar)]- Combining 
Eqs. (]4a|), (g) and (Q) we find the zero- momentum s-wave 
scattering length a 65 , which is plotted in Fig. H, together 



The ABC approximation is clearly much better than 
the DIS approximation. To obtain the DIS result, the hy- 
perfine level splitting is negelected, which is equivalent 
to setting k^ = in Eq. ([la|), resulting in Re(a 6 s) = 
aT cos 2 9 + + as sin 2 9+ and Im(a6s) = 0. For hydrogen 
and deuterium, terms containing k^}p,(R — ax,s) can in- 
deed be neglected compared to 1 at lower magnetic fields, 
so that DIS is a rather good approximation in that case. 
However, for 6 Li with its anomalously large value of aT, 
this is not the case, and one must then take terms of or- 
der fcgj (R — a^ into account. Also, for heavier alkalies, 
the DIS approximation becomes increasingly bad even at 
lower magnetic fields, due to the dependence of k^/g, on 
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the atomic mass and hyperfine constant. 

It is straightforward to show that at large magnetic 
fields, the ABC approximation gives Re(a6s) = ay, and 
Im(<265) = 0, independent of the choice of R. At B = 0, 
we find Re(a@5) = 3as — 2R ~ 57ao for R = 40ao, and 
Im(<i65) = — l/(k^ tan 2 9+) ~ — lOOao, in surprisingly 
good agreement with the results of the exact CC calcula- 
tion. Note that the imaginary part becomes rather large 
and contributes significantly to the elastic cross-section 
47r|a65 1 2 a t lower magnetic field. 



tively easily and accurately with the asymptotic bound- 
ary condition approximation. 
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FIG. 4. Exchange decay rate G65^6i of the Li gas. The 
solid line gives the CC result. The dashed lines represent 
1) the DIS approximation, 2) the ABC approximation with 
R = 40ao, and 3) the ABC approximation with R = 60ao. 

Similarl y, o ne can find the decay rate constant Gg5^6i 
from Eq. ( |4b|) combined with Eqs. (||) and (§). The re- 
sults are plotted in Fig. [|. Again, agreement with the 
exact CC calculation is very good, whereas the DIS ap- 
proximation is in error by several orders of magnitude. 
We believe that application of the ABC approximation 
to other alkali atoms for which the potentials are not so 
well known, can also give an accurate estimate for the 
collisional properties of interest pl|. Of course, the re- 
sults will vary slightly with the actual choice of R. It is 
difficult to predict a priori which R reproduces the ex- 
act results for an arbitrary atom best, but any R chosen 
as indicated previously, will give at least a good order 
of magnitude estimate. The DIS approximation on the 
other hand, can be off by many orders of magnitude. 

In summary, we have calculated the exchange decay 
rate constants and s-wave scattering length for four dif- 
ferent combinations of hyperfine states of Li by means 
of a coupled channels calculation. The scattering lengths 
attain their large and negative value of — 2160a only for 
B > 0.05 T. We found that the |{21}) combination has 
Feshbach resonances at B = 0.08 T and B — 1.98 T. In a 
magnetic trap, the lifetime of the |{64}) combination is 
slightly higher than the one for the |{65}) combination. 
The results of the exact calculation can be obtained rela- 
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